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We argue that supersymmetric flat direction vacuum expectation values can decay non- 
perturbatively via preheating. Considering a toy U{1) gauge theory, we explicitly calculate the 
scalar potential, in the unitary gauge, for excitations around several flat directions. We show that 
the mass matrix for the excitations has non-diagonal entries which vary with the phase of the flat 
direction vacuum expectation value. Furthermore, this mass matrix has zero eigenvalues whose 
eigenstates change with time. We show that these light degrees of freedom are produced copiously 
in the non-perturbative decay of the flat direction vacuum expectation value. 
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I. INTRODUCTION 

The scalar potential of the Minimal Supersymmetric 
Standard Model (MSSM) possesses a large number of 
F- and D-flat directions along which the scalar potential 
nearly vanishes [l], ^ . These flat directions can have im- 
portant cosmological consequences, including the genera- 
tion of the baryon asymmetry of the Universe through the 
out-of-equilibrium CP violating decay of coherent field 
oscillations along the flat directions themselves [1, |3, @ • 

Recently, much interest has focused on the cosmo- 
logical fate of flat direction vacuum expectation values 
(vev)s. In particular, it has been argued @ that in real- 
istic supersymmetric models, large flat direction vevs can 
persist long enough to delay thermalization after inflation 
and therefore lead to low reheat temperatures. Further- 
more, it has also been asserted that large flat direction 
vevs can prevent non-perturbative parametric resonant 
decay (preheating) of the inflaton since the inflaton decay 
products become sufficiently massive preventing preheat- 
ing from ever becoming efficient. These arguments hold 
so long as the flat direction vevs do not rapidly decay 
- they must persist long enough so that they can delay 
thermalization and block inflaton preheating. In [s] it 
was claimed that non-perturbative decay can lead to a 
rapid depletion of the flat direction condensate and thus 
precludes the delay of thermalization after inflation. It 
was also concluded that in order for the flat direction to 
decay non-perturbatively the system requires more than 
one flat direction Q. Finally, in Q it was pointed 
out that even in the presence of multiple flat directions, 
some degree of fine-tuning was necessary to achieve flat 
direction decay. 

An important aspect of this discussion centers on the 
issue of Nambu-Goldstone (NG) bosons. In general, su- 
persymmetric flat directions are charged under the gauge 
group of the MSSM. Consequently, the flat direction vev 
will break some or all of the gauge symmetries of the the- 
ory and thus we expect the presence of the associated NG 
bosons. In calculating non-perturbative flat direction de- 
cays, [8j considers a gauged U{1) model and constructs 
the mixing matrix for the excitations around the flat di- 
rection vev. The results in ^] show that in the single flat 
direction case, non-perturbative decay proceeds solely via 
a massless NG mode as only the NG mode mixes with the 



Higgs and all other massless moduli remain decoupled. 
Since the NG boson represents an unphysical gauge de- 
gree of freedom, it was concluded that no preheating 
occurs in the single flat direction case. As the appearance 
of a massless NG boson in the spectrum is a gauge depen- 
dent artifact, it remains unclear if the conclusions drawn 
about the system hold in the unitary gauge. In order to 
determine if flat direction vevs decay non-perturbatively 
into scalar degrees of freedom, the effect of the NG bo- 
son mixing in the scalar potential must first be removed. 
The process of removing the NG modes by switching to 
the unitary gauge changes the form of the mixing matrix 
among the left over scalar degrees of freedom. 

In this letter we consider toy models to demonstrate 
that, in the unitary gauge, the mixing matrix of the ex- 
citations around a flat direction vev permits preheating. 
Moreover, we find that flat direction decay depends on 
the number of dynamical, physical phases appearing in 
the flat direction vev. Specifically, a physical phase dif- 
ference between two of the individual field vevs making 
up the flat direction is needed. 

The outline of the rest of this letter proceeds as fol- 
lows: first we explicitly construct - in the unitary gauge 
- the mass squared matrix arising from the D-terms of 
a toy gauged U{1) model with three charged chiral su- 
perfields. We then present the formalism of preheating 
with multi-component fields and show that preheating 
occurs for the light moduli associated with the flat di- 
rection. We then analyze the specific dynamics of the 
background field equations for the toy models examined. 
Finally we evolve the background field equations for one 
of the toy models to obtain quantitative results. 



II. TOY MODEL WITH A GAUGED U{1) 
SYMMETRY. 



As an example, we examine a toy model which demon- 
strates the important features of supersymmetric flat di- 
rection vev decay. We introduce three complex super- 
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fields^ , $1, $2 and $3 charged under a U (1) gauge group 
with charges qi — +1/2, q2 = —1 and (73 = +1/2 respec- 
tively. The Lagrangian reads 



with 



1 



1=1 



(1) 



where = 9^ — iqiA^^ denotes the covariant derivative. 
The potential we consider arises from the supersymmetric 
D-terms and has the form 

^= y(9i|'i>i|' + 92|$2|' + <Z3|<i>3|')', (2) 

where g is the gauge coupling of the U{1) gauge sym- 
metry. In the above, we have neglected contributions 
from supersymmetry (SUSY) breaking and from any 
non-renormalisable terms arising from the superpoten- 
tial. These contributions are highly model-dependent 
and cloud the analysis we wish to present. A fully re- 
alistic model must include these additional contributions 
which can significantly affect the resulting particle pro- 
duction. The effects we investigate here do not depend 
on their inclusion in the quadratic part of the potential 
and so for clarity we neglect thcm^. 

The potential in eql2] admits several flat directions. 
Choosing one particular direction and including excita- 
tions around the vev we can write 



$1 
$2 
$3 



(^ + ei)e^(^^+v), 
(^ + e2)e'("^+^\ 



(3) 



where CTi 2,3 represent time dependent phases of the vevs, 
(p denotes the vev's time dependent amplitude^ while 
^1,2,3 and ^1,2.3 parameterize the six real scalar degrees 
of freedom corresponding to the excitations around the 
vevs. Note that the flat direction vev breaks the C/(l) 
gauge symmetry. Thus, out of the six real scalar degrees 
of freedom we expect one massive Higgs field and one 
massless NG boson, leaving four massless scalar degrees 
of freedom. 

The kinetic terms for the scalar fields play an impor- 
tant role in this analysis. Their expansion in eq[T]includes 
the term 



C D -ip^Aoidi ~ 2(72 + fT3) 



(4) 



which has the form of a coupling between the gauge 
field and the background condensate. Terms of this type 
will feed into the equations of motion for the gauge field 
which, in turn, will have an effect on the equations of mo- 
tion for the scalar excitations. By making a U{1) gauge 
transformation on the vev of the form 



(5) 



^ We use $ to denote both superfields and scalar components of 
superfields. 

^ SUSY breaking terms are needed for the computation of the evo- 
lution of the flat direction vev. 
^ Throughout this analysis we assume ip <C I/Sct. 



A 



2(72 - 0-1 - 



(6) 



we can gauge this term away and avoid a complicated 
analysis of the kinetic terms. The resulting form of the 
vev reads 



(7) 



where 7 = (cti — o'3)/2 and a = {ai +<J2 + o'3)/5 represent 
the two remaining independent physical phases. Follow- 
ing Kibble jT3j |. we can write the fields in the unitary 
gauge as 



$1) 




$2) 




3'3> 





^1 = ((^ + a)e'^"+''+^^+^\ 
*3 = {f + ^3)e 



(8) 



where ^1,2,3, and 9' denote the physical excitations - 
the NG boson has been removed^. We choose the par- 
ticular combination of field excitations appearing in eq[5] 
(the exponent in particular) in order to retain canonically 
normalised kinetic terms. 

On substituting the fields of eqlS] into the La- 
grangian given in eqlT] and defining the vector S = 
{£,i,£,2i£,3,(),()')'^ , we find the quadratic terms 



Ie^M'E+^-E^UE + ... (9) 



where the ellipses denote higher order terms and inter- 
actions. The matrix U given in the last term in eqlS] 
reads 



U 



/ 











<T + 7 


(T + 7 

V3 


\ 
















a 
.V3 















— (T+7 

V2 


(J — 7 
~73 






o"+7 
V2 





(T — 7 

-cr+7 
V2 










V 


cr+7 
V3 


a 

V3 








/ 



(10) 



while the mass matrix for the physical excitations ap- 
pears as 



/I -2 1 0\ 

' -2 4 -2 * 
1-2100 



V 0/ 



= BMlB^ (11) 



M| = Ml = 



with eigenvalues = 6(3(^9)^, M| 
Af| = (the entries of the diagonal matrix J^d)- B is an 
orthogonal matrix which diagonalizes A4'^ and Mi corre- 
sponds to the mass of the physical Higgs field associated 



* This can be verified by expanding out the scalar kinetic terms 
which reveals the absence of terms of the form 9'' (...). 
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with the spontaneous breaking of the U{1) symmetry. 
The four zero eigenvalues correspond to the massless ex- 
citations around the flat direction vev. 

The last term in eql9] appears as a consequence of 
the time-dependence of the background - it represents 
a mixing between the fields $1,2,3, d, 6' and their time- 
derivatives. The effect of these terms on the system be- 
comes clear if we make field redefinitions that remove 
the mixed derivative terms. The resulting transforma- 
tion leaves the system in an inertial frame in field space 
and leads to a time-dependent mass matrix. Defining 
S' = AE. {A is orthogonal), we find the condition that A 
must satisfy, in order for all the mixed derivative terms 
to cancel, to be 



A^A = u. 



The Lagrangian for the S' system now reads 



1, 



"/1 2 



1 



(12) 



(13) 



where TW'^ = AM^A'^ = ABMjB'^A^ = CMlC^ and 
C = AB. The matrix C is an orthogonal time-dependent 
matrix, with columns corresponding to the eigenvectors 
of A^'^. We now have a system of scalar fields with 
canonically normalized kinetic terms and time dependent 
eigenvectors. 

The central point of this discussion centers precisely 
on the appearance of the time dependent eigenvectors for 
the five scalar fields. This satisfies a necessary but not 
sufficient condition for preheating. In the next sections, 
we investigate the details of the non-perturbative pro- 
duction of the light scalar fields following the analysis of 



III. NON-PERTURBATIVE PRODUCTION OF 
PARTICLES 



Including gravity, the dynamics of the re-scaled con- 
formally coupled scalar fields, Xi = ^■=■^1 where a denotes 
the scale factor and E[ the z-th component of the vec- 
tor S' , are governed by the following equations of motion 
(sum over repeated indices is implied). 



(14) 



where dots represent derivatives with respect to confor- 
mal time t, and 



(15) 



where k labels the comoving momentum. Using an or- 
thogonal time-dependent matrix C(i), we can diagonalize 
fiy via C"^(i)r2^(i)C(<) = '^^(^), giving the diagonal en- 
tries (^j{t). Terms of the form ~ ^iJX arising from the 
kinetic terms do not affect the evolution of the non-zero 
k quantum modes . 

Once we have identified the basis in which the Hamilto- 
nian appears diagonal (via the orthogonal matrix C(t)), 
the study of particle creation by the time-varying back- 
ground proceeds as in [l^, [l^ [l3|, which extends the 



results of [I^. Following [I^, we assume that flij ini- 
tially evolves adiabatically by assuming that the ini- 
tial angular motion of the flat-direction varies slowly. 
This assumption allows us to define adiabatically evolv- 
ing mode functions with positive and negative frequency. 
We rewrite the quantum fields as mode expansions in 
terms of the mode functions and their associated cre- 
ation/annihilation operators which allows us to define 
the initial vacuum. During the evolution, the entries of 
flij do not necessarily change adiabatically and conse- 
quently we must find new mode functions that satisfy 
eqUD A new set of creation/annihilation operators re- 
quired to define the new vacuum can be related to the 
initial set using a Bogolyubov transformation with Bo- 
golyubov coefficients a and /3 (which denote matrices in 
the multi-field case). 

Initially a — I and /3 — while the coupled differential 
equations (matrix multiplication implied) 



a 



'iijja + ——(i — la — J(3 
2lo 



$ — — a + iuofi — Ja — If} 
2uj 



(16) 



govern the system's time evolution with the matrices I 
and J given by 

/ ^ \ (V^C^C^ + -1= C^CV^) , (17) 



J=-,[V^C-C^-^C-CV^y (18) 

Similarly to the single-field case it can be shown that 
at any generic time the occupation number of the ith 
bosonic eigenstate reads 

n,{t) = {I3*l3^)u. (19) 

As pointed out in [1, , there exists two sources of non- 
adiabaticity in the multi-field scenario. The first source 
arises from the individual frequency time dependence and 
appears as the only source of non-adiabaticity in the sin- 
gle field case. The second source appears from the time 
dependence of the frequency matrix giving rise to 
terms in ea ll6l proportional to / and J. This second 
source provides the most important contribution in our 
analysis and gives rise to non-perturbative particle pro- 
duction. 

Since initially a = I and /3 = 0, eqlH] shows that a 
non- vanishing matrix J is a necessary condition to obtain 
/3 ^ and hence ni{t) ^ 0. In general, we have 

C^C- = b'^A^AB = -B^UB (20) 

where A, B and U were defined in the previous section. 
For the toy U{1) example outlined above, J is a 5 x 5 
matrix in the Xi basis with non- vanishing components 



Ji,2 — J2.1 — 



2V3fc(fc2 -^M2)i/4 



7, 



(21) 



where Mi denotes the mass of the heavy Higgs field. 
These entries in the matrix J link the eigenstate of the 
Higgs {i = 1) with one of the light eigenstates {i = 2). 
We see that in the toy U{1) model, preheating can occur 
provided that 7 7^ 0. We address this point in a later 
section. 
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IV. MULTIPLE VEV AMPLITUDES 

We can extend the analysis of the previous sections 
by allowing the magnitudes of the individual field vevs 
to differ from one another. As above, we consider the 
case with three complex superfields charged under a J7(l) 
gauge group with charges qi = +1/2, q2 = —1 and qs = 
+1/2 respectively, and with the scalar potential given in 
eqlH We can write the flat direction with the following 
vev 



('i>i) 
('i>2) 



i=(^? + ^^)V^e-, 



(22) 



By substituting the above into the potential given in eqlH 
it can readily be shown that the configuration satisfies D- 
flatness. Expanding around this vev we have 



i (T1 + - 



$1 = {^1 + Ci)e 

<i>2 = [^-^{^l + vlf" 



6 e 



CT2 + 



(23) 



The resulting spectrum consists of one Higgs field with 
mass Ml — ?>g^{ip\ + ^2)1 ^'^d four massless scalar fields. 

We proceed, as before, by diagonalizing the kinetic 
terms and evaluating the J matrix given in ea llSI The 
non-vanishing entries of the J-matrix are 



Ji,2 — J2,i — 



VP 



Ml 



2V3fc(fc2 + M( 



U/4 



(26) 



which demonstrates that in this case preheating can take 
place provided that 7 0. 

It is instructive to compare the two cases considered 
thus far. The first flat direction contained a single vev 
amplitude, the second contained two independent vev 
amplitudes. The final result, however, is the same for 
both cases. This demonstrates a simple property of flat 
direction vev decay: the determining factor is not the 
number oj flat directions present in the system, but the 
number of fields that have vevs. In particular, a neces- 
sary (but not sufhcient) condition for non-perturbative 
production of particles is the existence of at least one 
relative physical and dynamical phase between the field 
vevs that constitute the flat direction. 



DYNAMICS OF THE VEV PHASES 



where the fields ^1,2,3 and 6*1^2,3 represent the excitations 
around the vevs. As in the previous case, we can use 
a gauge transformation to remove a phase from the vev 
structure that ensures the absence of terms of the form 
appearing in eqUl The form of the vev in this case be- 
comes. 



($1 



1 



2^,l/2gio- 



($2) = -^((p? + ^^) 

($3) = ip2e 



where a and 7 represent two independent phases^. 

In the unitary gauge, a form that preserves the canon- 
ically normalized kinetic terms reads. 



$1 = 



((pi+a)e V '^'^ 



(i=(^? + ^^)V2 + 6) 

((^2 +6)e 



e'(2/3)^/^ 

(^f+^i)i/V,(25) 



where ^1.2,3, S and 9' label the physical excitations 
around the vev once the NG boson has been gauged away. 



We now demonstrate that both physical phases are 
in general dynamical. In our particular toy example 
the cancellation of [/(l)"^ and mixed J7(l)-gravitational 
anomalies requires that we extend the field content of 
our model by including three additional complex super- 
fields, $4, $5 and <I>6. We assign the C/(l) charges and 
i?-Parity {Rp) as follows: 







^1 ^2 


$3 


$4 ^5 


$6 


(24) 


U{1) 


1/2 -1 


1/2 


-1/2 1 


-1/2 




Rp 


+ - 


+ 


- + 





This choice of Rp assignments forbids the superpoten- 
tial term $ii>2$3, thus preserving F-flatness. There ex- 
ist several possible flat directions for this particular field 
content. We assume that vevs for only $1, $2 and $3 are 
turned on, leaving $4, $5 and $6 with no vevs. With 
this assumption, the lowest dimension gauge invariant 
operators which have vevs are 



Oi = $i$2$3 



O2 



$f$2 



and O3 



<f2n. (27) 



Note that the last two operators depend on both phys- 
ical phases, a and 7. Using these operators, the lowest 
dimension terms appearing in the scalar potential, which 
are Rp invariant and phase-dependent, arise as soft SUSY 
breaking A-terms and appear as 



Again we have applied the limit <^ tpa, 1/37. If we do not 
apply this limit the gauge transformation parameter (A) needed 
to remove the linear term in Aq can be found by integrating 
the coefficient of the Aq term with respect to time. This is in 
general complicated and we choose to assume that ip is varying 
very slowly with time. 



V D 



E 



A, 



M3 



(28) 



where M denotes the cut off scale of the theory (e.g. the 
Planck mass or GUT scale) , nig represents the scale of the 
SUSY breaking, and Aij label dimensionless coefficients 
of order one. Lower order phase-independent interactions 
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will also contribute to the lifting of the flat direction and 
have the generic forms 



hi 



(29) 



where mf denote the soft SUSY breaking masses, and 
the second terms arise from loop corrections with Xij ~ 
g'^ml/ip'^ (see for example [s^ for similar loop induced 
terms). The potential for the single flat direction am- 
plitude case considered in section [TTl using the vev form 
shown in eq[71 becomes 



V D 



3 , „2 



A' 4 



2 " 4 

where A'^^ and A' denote combinations of the couplings 
discussed above. The potential for the multiple vev am- 
plitude case will be very similar with the obvious changes 
of vev amplitudes. The phase-dependent terms in the 
potential provide non-trivial dynamics for the phases a 
and 7 and will in general lead to 7 ^ and therefore 
a non-vanishing J matrix. As discussed above, the ap- 
pearance of a non-vanishing J matrix can lead to the 
non-perturbative production of particles by the rotating 
flat direction: the condensate can decay via preheating. 



VI. TWO INDEPENDENT FLAT DIRECTIONS 

A further instructive toy model consists of two inde- 
pendent flat directions existing simultaneously. We con- 
sider four chiral superfields ^1,^2,^3, and $4 charged 
under a gauged U{1) symmetry with charges ±qi and 
±52 respectively. The potential arising from the D-terms 
reads 



V^^iqi |$i|' - qi |$2|' + 92 l^sl' - 92 |$4|')'. (31) 



Although this toy model has been examined previously 
in Q , applying the methods outlined in the first sections 
of this letter helps establish the important properties of 
the model. The potential in eqlST] admits flat direction 
vevs of the following forms 



($1) = v'le'^S 
($2) = Vie'"% 



(32) 



($4) 



(p2e 
ip2e 



As before we can make a gauge transformation and re- 
move one phase in such a way that terms of the form 
shown in eq|4] vanish. The final form appears as 



j(cti+7-£2-) 

i(o-2+7-^) 

(y52e ' <!2V2 ' J 

($4) = <^2e'^"'"^^'. 



('fi) 

(*3) 



(34) 



demonstrating the existence of three physical phases. 
Transforming in to the unitary gauge, we can write the 
excitations around the vevs as 



2 = (¥^1+6)6 ^"1^1^ v^^l 



$3 = iV2+^3)e^^''^^^''^^ 



I 92 2 \ 

(35) 



$4 = (^2 + ii)e 



i(<y2+i- 



r) 



where tp' = \/2{qiLp\ + g2'/'2)^^^- The spectrum in this 
case consists of one massive Higgs particle and six mass- 
less scalar fields (the NG has been gauged away). Again, 
we must diagonalize the kinetic terms. Applying the nec- 
essary field redefinitions we are able to evaluate the J 
matrix. The non-vanishing J matrix elements read 



Ji,2 — J2.1 



^^1,3 — J3.I — 



Jl,4 — Ji,! — 



k-^WTMl i^2 



{CTI - a2)<f\ip2 



k - ^Jk'^ + Ml 7(^1 



(36) 



2fc(fc2+M2)l/4 ip' 



which depend on the relative phases between the field 
vevs. We should point out that only the Higgs eigenstate 
{i = 1) is distinguishable. The other indices label the 
light fields which at this level are all massless. Preheating 
is again possible provided two of the phases have non-zero 
time derivatives. Using the particular case with qi = 
92, we can write the scalar potential (see Appendix for 
details) yielding the terms 



We can write the excitations around the vevs as 



8 8 M^i^2 

Aims 22 2i(rTi+<T2-7^^^^) , /o'7^ 

Tm""^^^^^ (37) 

Clearly, non-trivial dynamics exist for the phases 7, ci 
and U2- 



$1 - (^i+ei)e^^''^+^\ 

$2 = (^i+6)e^^*^+-\ 

$3 = {V2+^3)e^''^^^\ 

$4 = (^2+^4)6*^''^ + ^^ 



(33) 



VII. NUMERICAL ANALYSIS 

As a proof-of-principle that achieves quantitative re- 
sults, we numerically analyse the model described in sec- 
tion IVII We use a simplified version of the potential 
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appearing in eql371 confining ourselves to the potential 



V 



>V2 



Aim, 4 ,4,^ 
8 

8 



h.c (38) 
This poten- 



where m^^ = m{ + mj, = + 
tial decouples the equations of motion for 7, cti and (72 ■ 
The equation of motion for 7 reduce to 7 = 0, and with 
the choice of initial conditions, 7 = 0, the effects of 7 
on preheating are removed. Our simplified potential al- 
lows us to numerically evolve the classical evolution of 
the flat direction vevs and analyze particle production in 
a self-consistent background. We also make the simplify- 
ing assumption setting Ai = A2 = X-^ in eq[271 Again, 
we stress that we use this grossly simplified potential sim- 
ply to demonstrate the quantitative behaviour of the toy 
model class. 

Measuring the conformal time in units of t ^ ft with 
/ = gipig and using the re-scaled flat-direction vev am- 
plitudes 




4 6 8 10 12 14 
Dimensionless conformal time 



FIG. 1: Occupation number for one of the excited fields as 
a function of dimensionless conformal time, obtained using 
ea ll9l after numerically integrating the background field equa- 
tions and Bogolyubov matrices; k = /ii/3 x 10"", other param- 
eters as explained in the text. The solid lines represents pre- 
heated fields with Lp2g /vio — 1 while the dashed lines indicate 
the preheated fields with 9220 /v'lo — 0-1 



a 



^2 



<^2o 



(39) 



we find the background equations. 



2 9 

Mia 



'2 



Fi + —^cos{4ai)^0 



F'' 



9 9 



Fo 



XFi 



^ ) cos (4ct2) = 



(40) 



where a prime represents a derivative with respect to r 
and 

a'; + 2a[^~^F^smiAa,)^0 



Fi 2g 



a'^ + 2a',-^-—[r^ Fi sin(4a2) = (41) 

describes the motion of the flat direction vevs; /ii = 
m^p-^/ f, ^2 = ^ip-il ! ■ The scale factor evolves as, 



t'2 



rAi^iFt 



a2 



cos(4cri) 



M2 



F' 



X 



A 



F^ 

— ^ C0s(4cr2) 



Initially the flat direction vevs correspond to a con- 
densate of coherent particles with vanishing momentum. 
The motion of these vevs, described by eqHOl and cgHTl 
and the interactions described in the previous section, 
cause the rapid decay of this condensate into a decoher- 
ent state of particles. FIG. 1 shows the occupation num- 
bers, ni{t), of these light particles as a function of con- 
formal time: the exponential growth of these functions 
signals the exponentially fast decay of the flat-direction 
vev. The two line types, solid and dashed, represent the 
ratios ^2qI^\o — li 0-1 respectively. We see that preheat- 
ing occurs over a wide range of the ratio (^20 /'/'lo- this 
numerical example, preheating effects do not vanish until 
^2qI^\o ^5 10^^. FIG 2 displays the resulting spectrum 
for one of the light fields, we see that production of higher 
momentum modes becomes kinematically suppressed. 



(42) 



where is the energy density of the inflaton field and 
/p = ifi^/Adpi is set to fp = 0.1 in our numerics. We also 
take /ii = 10^2, /i2 = 10^^/2, and X = p\ for computa- 
tional ease. As initial conditions, we start the flat direc- 
tion at rest, such that (</3i,2 exp(«(Ti^2))' — 0. We choose 
to set initially Fi^2 = 1, o'i.2 = 0.05, (t[ 2 — 0, a — 1 
and a' — fii, which implies Fl 2 — o-' = fJ-i ■ While these 
initial conditions do not present a realistic case (where 
/i ^ 10^^^ and F{ 2 ^ Mi); they do provide a numerical 
proof-of-principle similar to 0]. 




FIG. 2: Occupation number as a function of comoving mo- 
mentum obtained as for FIG. 1, with (^2o/<^io = 0-5 at time 
t = 20. 

We must stress that effects of SUSY breaking terms in 
the Lagrangian eq[T] will significantly affect the amount 
of particle production produced by the rotating conden- 
sate. A mass term for the light fields translates into a 
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momentum shift in eqllSI and this corresponds to a kine- 
matic suppression of the modes [8]. A reahstic model 
involving MSSM flat directions will in general contain 
many SUSY breaking terms and non-rcnormalisable op- 
erators, creating large model dependencies in the precise 
determination of the momentum shift. Consequently, we 
leave such a study to future work. 

VIII. CONCLUSIONS 

The cosmological fate of flat directions provides a ma- 
jor ingredient for the history of the early Universe. Flat 
directions can provide mechanisms for generating the 
baryon asymmetry of the Universe and can play an im- 
portant role in reheating after inflation. Our analysis 
stresses the use of the unitary gauge in which the physical 
content of the theory becomes manifest. By transforming 
to the unitary gauge, complications arising from mass- 
less NG modes in the mixing of the excitations around 
the flat direction vev are removed. The mixing matrix 
in this gauge defines the mass eigenstates of the physical 
scalar fields and determines if non-perturbative decay is 
possible. Since the mass matrix in the unitary gauge can 
contain time dependent mixing among all fields, one of 
the necessary conditions for preheating can be satisfied. 

Two further crucial conditions for preheating in our 
analysis center on the existence of physical relative phases 
between the field vevs that make up the flat direction(s) 
and that these phases possess non-trivial dynamics dur- 
ing the early universe. The first of these conditions gener- 
ally becomes satisfied if the difference between the num- 
ber fields that acquire a vev and the number of broken 
diagonal generators is larger than one - every diagonal 
generator removes one unphysical phase. The second 
condition generally becomes satisfied if terms which ex- 
plicitly depend on the phase differences appear in the 
scalar potential. The existence of gauge invariant prod- 
ucts of background fields exhibiting this phase depen- 
dence represents the crucial ingredient and determines 
the phase dependence in the scalar potential. Once these 
conditions are satisfied, the flat direction condensate can 
decay non-perturbatively via preheating. 



IX. APPENDIX 

In this appendix we specify a simple model to justify 
the form of the potential given in eql371 We have the 
following field content with U{1) charges and Rp assign- 
ments, 





$1 $2 $3 $4 


U{1) 
Rp 


1-11-1 

- + - + 



The lowest dimension gauge invariant operators are 

Oi =$i$2, 0-2 =$3*4, 

O3 =$1$4, C4 =$2«'3- (43) 

The lowest dimension terms which are Rp invariant and 
phase-dependent arise as soft SUSY breaking A-terms 
and appear in the scalar potential as 

^ 2 ' 8 M 1 8 M 2 

i=l 

^Mra^Q. M^2±o\^... (44) 
8 M 3 8 M * ^ ' 

where the ellipses stand for other terms of the same order 
with different products of the gauge invariant operators, 
loop induced contributions and higher order terms in the 
(1/Af) expansion. Substituting the vevs given in ealMl 
we generate the potential terms given in eg 1371 
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